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ABSTRACT 



Composite strength is an inverse function of the size of the com- 
posite. As the use of composites expands into larger applications, 
such as airplane wings, missile components, and ship superstructures, 
the ability to accurately predict composite performance for large 
applications has become more important. The composite failure pro- 
cess is sequential and initiates with early breaking of the weak fibers. 
Concentration of breakage sites accumulates and leads to ultimately 
catastrophic failure. Prestressing fibers prior to solidification of the 
matrix has been demonstrated to increase the reliability of the com- 
posite by minimizing the spatial concentration of the breaikage sights. 

This study concentrates on quantifying the level of preload and 
gauge length to optimize the prestress effect. Computer simulations of 
graphite bundle tests were used to validate the data analysis method- 
ologies applied to actual AS-4 graphite bundle tests. The actual 
experimental results are consistent with computer-simulated 
behaviors. 
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I, INTRODUCTION 



In the 1980s, the capability of military equipment has increased 
exponentially. Modern aircraft fly further and carry greater payloads, 
missiles have larger warheads and greater ranges, tanks have bigger 
guns and tougher armor, and ships have more capable but lighter 
weapon systems along with improved armor plating. Engineering 
improvements in the areas of electronics and composite materials 
have been on the leading edge of this technological revolution. 

As reported in the fourteenth edition of Ships and Aircraft of the 
U. S. Fleet, the AV-8B Harrier II VSTOL attack aircraft “provides twice 
the payload of the AV-8A with up to 9200 lbs of external stores.” The 
primary difference between the AV-8A and the AV-8B is the extensive 
use of composite materials in the wing and fuselage. This is just one 
example of the use of composite materials in existing military applica- 
tions. New applications include the VS-22 Osprey tilt rotor aircraft, 
whose wing and fuselage will be entirely made of graphite epoxy com- 
posite material. Emerging “stealth” technology and its application to 
aircraft is a further new application of composites. 

Composites offer designers many attractive features. These 
include very high strength-to-weight ratios, high stiffness-to-weight 
ratios, wear resistance, controllable heat expansion coefficients, cor- 
rosion resistance, fatigue resistance, low price, and ease of manufac- 
turing and repairability. The field is expanding rapidly and research is 
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bringing to the designer stronger, tougher, and stiffer materials every 
day. 

A hurdle to be overcome for composites, however, is the general 
lack of understanding of the physics and mechanisms at work in the 
composite which give these materials their unique properties. This 
lack of understanding has restrained designers’ ability to accurately 
predict composite strength. To compensate for this, large factors of 
safety are being used in the design of composite structures, which in 
turn adds weight and size and causes other related problems that 
reduce design efficiency. As the application of composite materials has 
expanded to larger and larger components, such as entire fuselages, 
rocket motor cases, and gun turrets, this problems has become 
magnified. 

Composite materials are two or more dissimilar materials that 
when properly combined form a new material whose properties 
exceed either of the individual constituents. For the designer, com- 
posite materials provide an ability not only to optimize his design but 
also to optimize the material to be used for the design. Material prop- 
erties of conventional materials that can be improved upon by using 
composites include stiffness, corrosion resistance, wear resistance, 
strength-to-weight ratio, and radar wave absorbability and reflectivity. 
Composites in use today include fiberglass, WEST (Wood Epoxy 
Saturation Technique), graphite-epoxy, Aramid-epoxy, and Boron- 
epoxy, Composite materials also include graphite-aluminum and hybrid 
composites of more than two constituents. 
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Commonly accepted types of composites are the fibrous compos- 
ites, which consist of fibers in a matrix; laminated composites, which 
consist of layers of various materials; and particulate composites, which 
are composed of particles in a matrix. The composite is fabricated by 
combing the reinforcement material into a generally ductile matrix 
such as epoxy. This lay-up is called a lamina. Several layers of lamina 
may be combined in prescribed geometrical orientations to create a 
laminate. The designer specifies his materials and then designs the 
laminate to meet his requirements. 

Until the early 1980s material properties of composites were 
predicted by a simple principle called the “rule of mixtures.” This 
approximation assumed that each constituent of the composite con- 
tributed to the overall composite property in direct proportion to its 
property weighted by its volume percentage in the composite. This 
linear weighting method was satisfactory for composite physical prop- 
erties, which are averaged combinations of the constituent properties 
such as density and stiffness. 

However, the “rule of mixture” is very inaccurate for those physi- 
cal properties which are governed by extreme values such as perme- 
ation and strength. For such properties, considerable testing of a 
specific composite is needed to achieve representative material prop- 
erty values suitable for engineering design and reliability assurance. 
Not only are large amounts of data required for each specific 
composite but the designer’s material selections become limited to 
only those composites which have been tested. This characterization 
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methodology could not keep up with the rapid pace of the 
introduction of new materials. 

Clearly, a characterization method was needed to be able to 
predict composite material properties, particularly strength, if the full 
advantage of composites was to be obtained. Much work has been 
completed in the area by researchers such as Phoenix, Rosen, Harlow, 
and Wu. The most current papers on reliability and strength predic- 
tions of composites may be found in References 1, 2, 3, and 4. The 
most significant part of their efforts is that they have identified that 
the strength of composites is a probabilistic manifestation of the sta- 
tistical strength of the fiber strength. The statistics of the fiber 
strength govern the number of failure sites in a composite. The proba- 
bilistic spatial distribution of the failure sites together with the failure 
mechanism govern the strength of a specific composite. 

Different probabilistic models have evolved to capture the essence 
of the failure mechanisms. The simplest model is the “weakest link 
model,” which envisions the entire composite as a long thin fiber; the 
entire chain (the composite) fails when the weakest link of the fiber 
fails. This model does not account for the redundancy provided by the 
matrix and as a result grossly underestimates the strength of the com- 
posite. The “equal load-sharing model” accounts for the redundancy 
provided by the matrix through definition of an effective length, an 
isolation parameter for the initial breaks. However, this model does 
not account for the stress concentration around the breaks and, as a 
result, over-predicts the composite strength. Finally, the “local 
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load-sharing model” accounts for both the matrix redundancy and the 
local stress concentration. These refinements lead to the most 
realistic predictions to date. 

There are two phenomena relating to the statistical strength of 
composites. The first is that the larger a composite becomes, the 
weaker it gets. The second phenomenon is that the larger a compos- 
ite, the less strength scatter among such large structures. This is due 
to the the lower weak tail of the strength distributions of the con- 
stituent fibers. This effect has been demonstrated in the laboratory 
and in the field. 

An extension of this theory then becomes that if the fibers are 
preloaded to break the weak sections of the fiber prior to their lay-up 
into the matrix, then the lower tail of the composite strength 
distribution is reduced and, in fact, becomes bounded by a value. This 
process is called prestressing. For a prescribed value, the strength and 
reliability of the composite may then be theoretically calculated for a 
given fiber and the magnitude of preload based on the modified fiber 
distribution. The variables to be determined to prestress a composite 
include the magnitude of the preload, the gauge length to which the 
load is applied, and the modified fiber parameters. 

The ability to accurately recover the fiber strength parameters for 
small (less than 2.5 cm) gauge lengths is very important if the pre- 
stress effect is to be optimized. The bundle testing done in the labora- 
tory yields data that includes not only the fiber data but also machine 
compliance, slack, and friction. For short gauge lengths, the 
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contributions of slack and friction are small but compliance becomes 
relatively large. The portion of displacement contributed to by the 
experimental system compliance in a bundle test needs to be removed 
from the data if the true fiber parameters are to be obtained. 

This study entailed the mathematical simulation modeling of 
strength prediction for a bundle of fibers with no matrix. This simula- 
tion was run on an IBM Personal Computer using Microsoft Fortran 
4.01 for source code and Lotus 1-2-3 for graphing. When the simula- 
tion program had been validated, a zero gauge length bundle test was 
completed to measure the INSTRON fiber bundle testing machine 
compliance. A set of procedures was written to optimize a curve fit of 
the compliance data to provide an explicit equation which models the 
INSTRON system compliance. 

When this procedure had been validated, bundle tests of varjdng 
lengths were tested at the Mechanics of Materials for Composites Lab- 
oratory at the Naval Postgraduate School. The displacement con- 
tributed to the test data by the INSTRON system compliance was 
removed from the total displacement. Under another thesis effort, this 
rectified data was analyzed in order to extract the fiber statistical 
strength distribution parameters. 

The purpose of this investigation was to validate the experimental 
procedures and data analysis methods with the use of computer simu- 
lated data. Data generated from Monte Carlo simulation using known 
parameters is then utilized to assist in experimental design and the 
associated prediction theory. Experimental design and data 
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interpretation methodology provide a proven foundation toward iden- 
tification of an optimum gauge length for prestressing in order to 
guarantee the reliability of composite structures. 
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II. BACKGROUND 



The tensile failure process of composite materials is a complex 
sequential combination of many different processes. These processes 
are controlled by the statistical strength of the fiber and the relative 
spatial clustering of the broken fiber sites within a composite. Given a 
probabilistic model which is capable of characterizing the sequential 
failure mechanisms, a formula can be derived to predict the 
probabilistic strength of the composite in terms of the statistical 
strength of the fiber. The parameters for such a probabilistic models 
can then be expressed in terms of other common measures of central 
tendencies (e.g., mean composite strength) and dispersions (e.g., 
strength scatter). 

A. COMPOSITE STRENGTH PREDICTION 

The prediction of the probabilistic strength of a composite as a 
function of fiber statistics was first introduced by B. W. Rosen in 
Reference 5, in the early sixties. He refined the simple weakest link 
model (which is essentially a series model) by accounting for the 
redundancy introduced by the presence of a matrix binder. The com- 
posite geometry he examined is shown below: 

• a single layer of fibers 

• evenly distributed throughout the matrix 

• no fibers are in contact with other fibers 

• the load is purely axial and in tension. 
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The strength of a single lamina of composite material is con- 
trolled by four factors: 

• The strength and modulus distributions of the fiber. 

• The strength and modulus distributions of the matrix. 

• The interaction between the fiber and the matrix with respect to 
load sharing and stress concentration. 

• The location distribution of relative strengths of the fibers with 
respect to the location distribution of relative strengths of the 
matrix. 

When a composite is pulled in tension, because the fiber generally 
has a much higher modulus as compared to the matrix (e.g., graphite 
108 Psi versus epoxy 10"^ psi), the fibers will carry most of the load. 
The load-carrying capability of the matrix can be neglected. 

As the load is increased, each fiber will carry an equal amount of 
load that is equal to the applied load divided by the number of fibers. 

Pf=Pt/n (1) 

The load on each fiber will be equal, but because the diameter of 
the fibers is not constant, the stress carried at each point in the fiber 
will be different. The stress will vary inversely to the diameter of the 
fiber. Some of the variations are regions of gradual thinning, others are 
indentations or nicks. These areas of reduced thickness or indenta- 
tions or nicks are called “flaws.” These flaws reduce the load-carrying 
capability of the fiber by reducing its thickness £md introducing stress 
concentration areas in the fiber. 
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A great deal of study has been done in identifying and quantifying 
the effect of these flaws on a microscopic scale. The work of Phoenix 
and Harlow proposed that on a macroscopic scale, if the flaws are dis- 
tributed randomly with respect to location on the fibers and normally 
with respect to the severity of the flaw, then the strength of a fiber 
can be modeled as a long chain of segments. The strength of the seg- 
ments may be modeled using the Weibull distribution. The model 
selection is expanded in Appendix A. The Weibull distribution is an 
appropriate model for the probability of failure of a chain with many 
links. 

F(x) = 1 - exp I-(x/px)“] (2) 

Where x = is the load applied to the fiber (the chain) 

a = is the shape parameter characterizing the variability 

Px = is the scale parameter characterizing central location 
of the failure strengths. 

The Weibull distribution for strength predicts the “size effect.” 
That is, if the parameters are known for one given length, the 
parameters for any other lengths of the same material can be com- 
puted. This is also expanded in Appendix A. 

As the load increases on the material, the stress on each segment 
of the fiber will increase as well. When the load reaches a point where 
the resulting stress is equal to the ultimate strength of a segment, that 
segment of fiber will break. This break will release a certain amount of 
strain energy into the composite. It will also cause an immediate local 
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load concentration that will be felt by the adjacent and surrounding 
fibers. 

The matrix, via its shear stress carrying capability, Avill attempt to 
absorb the released strain energy and distribute the local load concen- 
tration to the surrounding fibers. The actual values of the stress con- 
centration factors caused by the fiber break are not known, but shear 
leg analysis yields reasonable approximations. 

Once one fiber segment breaks, many events can occur. The first 
is that the matrix is strong enough to absorb the released strain 
energy and the surrounding fiber elements are capable of carrying the 
local load concentration. If this occurs, the material is stable and the 
load on each fiber will increase to 

Pf=Pt/(l-n) (3) 

Where n represents the number of broken fibers. 

Another possibility is that the matrix can not absorb the released 
strain energy. In this case the matrix will split as it fails in shear. The 
material will continue to tear until its radius is such that the stress 
concentration decreases to a point below the ultimate strength of a 
section of the matrix. The fibers in this region will then be subjected 
to additional stress as the matrix pulls away. If the segments of the 
matrix are strong enough to carry the additional load, the material will 
then again become stable. 

A further possibility is that the adjacent fiber segment may not be 
able to withstand the stress caused by the first fiber failure. In this 
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case, that fiber segment will break and cause additional strain energy 
to be released. This reaction is called auto clustering. The process will 
continue until it is interrupted by a fiber segment that is strong 
enough to carry the stress or the entire piece fails catastrophically. 

Some observations of the description above: The composite may 
fail due to either the local strength of the matrix or the local strength 
of the fiber or a combination of both. The likely initiation of the failure 
process will occur with a fiber break because the strength scatter or 
variance of the fiber is typically greater than that of the matrix giving 
rise to nucleation sites at already high stress regions. The chain of 
events that occur after that first break is dependent on the strength 
statistics of the adjoining matrix and fibers. If there are many strong 
segments surrounding a weak segment that fails, then additional 
material will not fail. If, however, this is not the case, then further 
damage occurs which may lead to final unstable catastrophic failure. 

From the previously described failure process, it can be observed 
that the strength of the composite is dependent on the nucleating 
weak segments and the probability of there being another weak seg- 
ment nearby. Therefore, the shape factor of the fiber strength distri- 
bution is the critical element that controls composite strength and 
reliability. Because of the mechanism of failure, as a composite 
becomes larger, the number of the population of fiber segments 
increases and the likelihood of the existence of very weak segments 
becomes larger. This causes the strength of the entire composite to 
decrease. 
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Another observation is that because the strength of the composite 
is controlled by the weak fiber segments, the strength may not be 
solved for explicitly because the actual value of the weak segments is 
probabilistic. The strength of adjoining segments is probabilistic, as 
well. 

Prestressing solves two critical problems. First, it sets the lower 
limit of the strength of the segments by eliminating the weak seg- 
ments by causing them to break prior to the application of matrix. 
Second, prestressing causes the weak segments to be broken without 
local load stress concentrations affecting the adjacent fibers and with- 
out the localized dynamic strain energy release because the weak seg- 
ments are broken before the matrix is solidified. This results in auto 
clustering being minimized. 

The critical implementation problem is then to find the "segment 
length” to ensure that the majority of the weak segments are elimi- 
nated and then determine a prestress level. The prestress level may 
be optimized to a level where the advantages of breaking the weak 
segments is balanced by the introduction of clusters to the matrix. 

B. PRESTRESSING 

G. J. Mills and associates, under an United States Air Force-spon- 
sored project, conducted prestress tests on Boron fibers in the early 
1970s. In his reports, documented in References 5, 6. and 7. he 
detailed how his work demonstrated the feasibility of the concept and 
provided data from early experiments that supported his expectations. 
His paper did not attempt to model or predict the effect of 
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prestressing. He did note the effect of var 3 dng prestress gauge length 
and that as the gauge length decreased, the strength of the fiber 
increased, but he did not attempt to quantify an optimum length. He 
described the strength increase of the composite purely in terms of 
removal of weak segment and reducing the tail-end scatter of the 
fiber. He did not attempt to predict the performance of the composite 
once the matrix had been set. He did report that the strength did 
increase. 

In 1986, Lt. David Bell and Professor E. Wu of the Naval Postgrad- 
uate School conducted prestress tests on AS-4 graphite/epoxy. The 
results are published in Reference 9. The report’s experimental data 
confirmed the phenomenon of strength improvement after pre- 
stressing. This was in concurrence with Mills’ laboratory data. Both 
papers accounted for the improvement by theorizing that the process 
of prestressing improved the strength of a composite because it 
reduced the strength scatter of the composite. The work done by Bell 
and Wu was only done with one gauge length, 10", and even with the 
problems encountered in handling the samples, the material did show 
a significant reduction in scatter and therefore an increase in the 
usable strength reliability. 

Implementation of these findings to a practical application 
requires a model that can adequately predict the effect of prestressing 
as a function of prestress level, gauge length, and fiber strength 
parameters. 
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C. TESTING METHODS 

There are two methods in which fiber strength statistics can be 
measured. One method is to prepare and perform a breaking strength 
test on one individual fibers. The second method is to form many 
fibers into a bundle of parallel fibers and then perform a controlled 
strain rate failure test. 

The first method is very time consuming and error prone. The 
fibers are very thin (in the order of 10 |im) and handling requires 
delicate treatment. Laboratory conditions must be carefully monitored 
as the small fibers are susceptible to damage caused by temperature, 
wind, abrasion, and other disturbances. In addition, each test only 
provides one single point of data. If failure probability data is required 
down to .001 failure rate, then at least 10,000 samples should be 
tested. 

Obviously, the second method is the preferred method. The bun- 
dle test, performed in our case on bundles of 1,000 filaments, will 
yield 1,000 data points with one test. The handling problems, while 
still present, are on a much more manageable level. The data analysis 
of the load vs. displacement data Avill yield the strength distribution 
parameters for the fiber that is equal to 1,000 individual filament 
tests. 

It is important then that we measure the true load and displace- 
ment for just the fibers. The displacement that the testing machine 
measures is the sum of the fiber displacement, machine compliance, 
electronic noise, slack in the individual fibers, and friction between 
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the fibers. Figure 1 demonstrates that test data is the combination of 
all of these factors. 




Figure 1 

Typical Bundle Test Data Graph 



The compliance is due to the mechanical construction of the 
testing machine. The machine operates with gears, clutches, and 
motors. Each of these has an initial start-up slip or backlash. The 
machine also has grips which hold the samples and the bases which 
will elongate under load. They are very stiff and the elongation is small 
when compared to the bundles, except when the bundle gauge length 
is short. This is illustrated in Figure 2, where the displacement is 
shown as percent strain. The compliance for the long sample has a 
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greater slope than the short sample because the slope of the compli- 
ance in displacement multiplied is by the gauge length for percent 
strain. 





GL = Long (> 25 cm) GL = Short (< 2.5 cm) 

Figure 2 

Effect of Compliance on Long-Gauge Length 
vs. Short-Gauge Length 

The compliance of the machine can be measured if the there is no 
contribution to the displacement from the fibers. This is achieved 
using a zero gauge length sample. A sample must be utilized to include 
the sample grip mechanism in the compliance. A curve-fitting 
approximation is then applied to the data to create a continuous func- 
tion of d in terms of P. This function may then be subtracted from the 
test data to eliminate the effect of compliance. This procedure and 
testing is expanded in Appendix C. 
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The effect of noise is sensed in the load cell. The sensitivity of the 
load cell is 1 out of 10,000. This is due to the limit of the analog-to- 
digital converter, which can only manage four bits of data. The amount 
of noise in the machine sensed by the load cell is very small when 
compared to the sensitivity of the load cell. As the load increase and 
the output changes in order of magnitude, the noise is not even 
noticeable. The noise in the machine is considered to be random 
electrical white noise due to interference in the data transmission 
lines. Observation of the data indicates that the noise is not very 
significant. 

A simulation program was written during this study that allowed 
the user to input random noise up to any level. This will allow simula- 
tions to be run and investigate the effect of noise on the fiber strength 
parameters recovered from the simulated test. 

The slack is due to the limitation of the sample preparation that 
all the fibers do not have exactly the same gauge length. This problem 
is found to be greater -with the longer samples as compared to the 
short samples. The slack can not be eliminated by current sample 
preparation techniques, but extreme care by the technician in 
preparing the sample can minimize the effect. 

Lt. Joseph Schmidt is examining the effects of slack on the calcu- 
lated fiber parameters and expects to publish the results of his study 
in September 1988, from the Naval Postgraduate School, Monterey. 

The effect of friction in the bundle test cannot not be definitively 
characterized. The presence of friction is manifested by large vertical 
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drops in the load displacement curve, perhaps due to the tangling of 
broken fibers with unbroken fibers. The friction also will cause a fibers 
to break more than once. The effects of friction do increase with 
longer gauge lengths and are reduced with shorter gauge lengths. 

With these factors that affect the data of the testing machine 
accounted for. the true fiber strength statistics may be recovered from 
the bundle tests. These statistics may then by used to optimize the 
prestress effect by calculating the effect prestressing will have in 
altering the fiber strength distribution. 
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III. PRESTRESSING MECHANISM 



The prestressing process is to preload the composite fibers prior 
to solidification of the matrix to cause the weak segments of the fibers 
to break. The load that was carried by the broken filaments is then 
shared equally among the remaining fibers. Equal load sharing is in 
effect because the pre-solidified matrix cannot support a shear load 
and cause stress concentration cells to form. Upon solidification of the 
matrix, however, local load sharing is in effect, thereby increasing the 
local redundancy. Because the weak segments of the fibers were 
broken prior to solidification of the matrix, the local load concentra- 
tion stress cells are not as severe as if the fibers had broken while in 
the matrix. The reduction in the severity of the stress concentration 
cells provides the reliability enhancement of the composite and is the 
motivation for prestressing. 

The parameters that affect the prestress effect are gauge length, 
prestress magnitude, and the fiber strength distribution statistics. The 
optimum gauge length to prestress the fibers would be the ineffective 
length. The ineffective length, as described by Rosen in Reference 5, 
is a function of both the tensile strength of the fiber and of the shear 
strength of the matrix. This makes the ineffective length unique for 
each composite, and because both strengths are random variables, the 
ineffective length for a particular location in the composite is random 
as well. 
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As the gauge length increases, the likelihood of breaking all of the 
weak segments in the fiber decreases. This is because using the equiv- 
alent load sharing model, it is only possible to cause one break in each 
fiber over the gauge length. The break occurs at the weakest flaw, 
which may occur before reaching the desired prestress level. As illus- 
trated in Figure 3, this allows other flaws weaker than the prestress 
level to remain unbroken. 



r INEFFECTIVE _ PRESTRESS 

L LENGTH GAUGE LENGTH 



GL 









Flaw 

Broken 



Flaw 

Not Broken 



M=1 M=3 M=6 

Figure 3 

Effect of Preloading for Various Gauge Lengths 
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When the prestressed fiber is placed into a matrix, the local load- 
sharing model becomes applicable. With this model, the weak flaws in 
the composite will break under load and generate the larger stress 
concentration cells which prestressing is designed to minimize. The 
mathematical modeling of solving for the probability density function 
of remaining weak flaws given gauge lengths greater than the ineffec- 
tive length and prestress levels is very complex. This problem is being 
treated in a separate project. 

The preload level has two effects on the prestress effect. For 
higher prestress levels, more fiber breaks are introduced and more 
stress concentration cells will be created in the the composite once 
the matrix has solidified. The severity of the concentration cells is a 
function of the ineffective length. The shorter the ineffective length is, 
the smaller the stress concentration cell is because of less local load 
sharing. Second, the higher the preload level, the more difficult it 
becomes to grip or maintain the desired gauge length. The fibers by 
their nature have a small diameter, on the order of 5 microns for 
graphite, and are very susceptible to handling damage. 

The underl 3 dng fiber statistics also have an effect on the prestress 
effect. The fiber strength variability (the Weibull model shape parame- 
ter a), has the greatest influence on the prestress effect. As seen in 
Figure 4 , if the shape parameter is very high (a > 20 ), i.e., the fiber 
has low variability, then the prestress effect will not be very helpful. 
This is because the range of the distribution of the fiber strengths is 
very small and the weak tail that prestressing eliminates or reduces is 
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inherently not present. Fibers of this strength uniformity, however, 
are not available at this time. If the shape parameter is small, then the 
scatter is much greater. With this scatter, the weak tail may be altered 
or removed by prestressing. 




Figure 4 

Strength Distribution for Fibers (Weibull Model) 

The quantification of the prestressing effect can be written as a 
function. 



PS = f (<D, 'F, a, p) (4) 

where O = Gauge Length 
'F = Preload Level 
a = Fiber shape parameter, Weibull 
P = Fiber scale parameter, Weibull 
The functional interrelation between these functions is presently 
not known and appropriate probabilistic modeling will be required. 
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One approach is to optimize each variable separately where possible. 
For example, the Weibull strength parameters characterize the fiber 
strength distribution. The prestress effect of the preload level and the 
gauge length will alter these parameters. The optimization of the 
effect on the alteration of these parameters on the reliability of the 
composite can be investigated as one means of optimizing the pre- 
stress effect. 

The probability of failure (F*) vs normalized load (x*/p) graph for 
the fibers (F) and composite (C) with the Failure Probability Density 
Function superimposed is shown in Figure 5. This figure demonstrates 
how the PDF can be transposed to a linear curve. The linear curve is 
much easier to use and because it reveals in great clarity what is hap- 
pening at the weak tail. This graph will be used to show the effect of 
prestressing with respect to failure probability. 




Probability of Failure vs. Normalized Strength 
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The probability of failure for the composite (C) is calculated from 
the fiber parameters. The composite reliability will diverge from the 
fiber reliability at Xc because of the local load-sharing stress-concen- 
tration cells. These cells cause the composite fibers in the matrix to 
break earlier because the fibers adjoining a broken fiber have to carry 
the load carried by broken fiber by themselves, not spread equally 
amongst all of the remaining fibers. Until the load Xc, the two curves 
do not diverge and the reliability of the composite is equal to the reli- 
ability of the fibers. 

If the fibers have been damaged so that their distribution curve is 
modified to appear like that shown in Figure 6, the result on the com- 
posite strength can be seen. As the composite graph swings upward, 
its reliability is decreasing. The composite is stronger than the fiber 
due to the load sharing but its reliability for loads greater than Xc has 
diminished. The composite reliability for loads below Xc has not 
changed at all. 

If the fibers’ strength distribution could be improved, the failure 
probability curve would be shifted, as shown in Figure 7. The impact 
on the composite strength is that the curve swings downward. This 
downward swing reflects greatly increased reliability. In fact, the 
composite will not fail theoretically at all for any load less than Xc. 

A method to achieve the fiber distribution in Figure 7 would be to 
proof test all of the composite fibers up to Xc and then remove any 



25 



O) 

3 



03 



O 

L- 

Cl 




Normalized strength 



Figure 6 

Probability of Failure vs. Normalized Strength (Damaged Fibers) 
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Figure 7 

Probability of Failure vs. Normalized Strength (Improved Fibers) 
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broken fibers from the bundle prior to adding the matrix. The effect of 
this process on the failure probability curve is shown in Figure 8 as 
curve R. The cost of this operation, considering the extreme thinness 
of the fibers, the susceptibility to damage due to handling, and the 
extreme amount of fibers that make up the composite, make this 
process in reality impossible. 



Xp - Proof test Fibers 

f “ Original distribution f 

R - Remove broken fibers 
Resultant distribution 

NRL - Not Remove broken fibers 
Long Gauge Length 
Resultant distribution 

NR5 - Not Remove broken fibers 
Short Gauge Length 
Resultant distribution 



Figure 8 

Effects on Failure Probability Curves Due to Proof Testing 

Instead, what may be done is to preload the fibers to Xc on a gauge 
length equal to the ineffective length of the fiber. Because of the han- 
dling difficulties, the broken fibers will not be removed. The effect of 
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this process on the failure probability curve is shown in Figure 8 as 
curve NRS. 

The NRS curve has swung downward as compared to the original 
fiber curve, but it has not swung as far down as if the broken fibers had 
been removed. The reason it does not swing as far is the local load- 
sharing stress-concentration cells that are created when the matrix is 
added. 

If the fibers are preloaded to Xc but on a gauge length that is 
greater than the ineffective length, the effect on the failure probability 
curve is shown on Figure 8 as curve NRL. This curve swings downward 
from the original fiber curve even less because of the probability that 
some fibers remain with segments that are weaker than the preload. 

As previously stated, the strength of the composite may be calcu- 
lated from the strength of the fibers. If the the fiber strength distribu- 
tions are altered as shown in Figure 8, the composite strength 
distribution will be altered in the same direction as well. Using this, a 
designer of a composite may then use a failure probability curve to 
determine whether his design meets his reliability specifications. If it 
does not, he may then decide to improve the composite by prestress- 
ing the fibers to a critical load to alter the failure probability to come 
within the requirements. 

The movement of the failure probability curve can be predicted by 
probability modeling based on the preload level, the gauge length, the 
ineffective length, and the original fiber strength parameters. The 
results of such an analysis (by Wu and Harlow) on the fiber swing for a 
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fiber preloaded at the level 0.3 * P , for an a of 5.0 and prestressed at 
multiples of its ineffective length, are shown in Figure 9. The graph 
demonstrates that even if the prestress gauge length is very large 
compared to the ineffective length, significant improvement of fiber 
strength distribution can be achieved. 

The optimization of the prestress effect is to identify the ineffec- 
tive length or to find the shortest prestress gauge length that can be 
applied. 
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Figure 9 

Probability of Failure vs. Normalized 
Strength Prestress Effect Simulation 
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IV. EXPERIMENTATION 



The experimental work was directed in two areas. The first was to 
measure the machine compliance of the INSTRON 4206 { a universal 
materials testing machine) in order to characterize the testing 
machine system compliance denoted by the function 6c (P). Once this 
was completed, the second portion of the experimentation was to 
perform failure tests on AS-4 Graphite bundle samples of various 
lengths in order to determine the strength characteristics and to 
observe the effects of gauge length on the test results with respect to 
compliance, slack, noise, and friction. Bundle samples were tested 
with no oil and then with oil in an attempt to reduce the effects of 
friction, particularly on the longer samples of 50 and 25 cm. Small- 
gauge lengths of 2.5 and 1 cm were tested with particular interest 
with respect to the examination of their suitability to preloading. 

A. COMPLIANCE TESTING 

The compliance testing was accomplished using a sample illus- 
trated in Figure 10. The copper tabs were prepared and washed with a 
dilute acid to enhance the performance of the adhesive. The graphite 
bundle was glued to the lower set of copper tabs with acetylene adhe- 
sive. The upper tabs were then set so that the distance (d) between 
the the upper and lower tabs was as near zero as possible. The 
graphite bundle was then glued into place with the same adhesive. The 
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sample was prepared by Composites Laboratory Technician. Jim 
Nageotte. 



Adhesive - 

Graphite 

Bundle 

Adhesive ■ 



Copper Tabs 




Copper Tabs 



Figure 10 

Compliance Test Sample 



The tensile test was performed on the INSTRON 4206 testing 
machine and associated INSTRON software, version 4.01 . A plot of the 
load vs. displacement is shown as Figure 11. The data and test sum- 
mary forms generated by the software are enclosed in Appendix E. 

The compliance data was then fit to the function shown as Equa- 
tion 4. This was done in accordance with the procedure for curve fit- 
ting detailed in Appendix C. The curve fitting was accomplished with 
the use of LOTUS 1-2-3 spreadsheet program and an IBM/AT. The 
resulting coefficients for the curve fit of the compliance are listed in 
Table I. 



32 



OCL 10508 ALL 




odd 



Figure 11 

Compliance Test Load vs. Displacement Graph 
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Displacorngnt (mm) 



INSTRON 4206 Compliance Curve Fit Function 



ForO<5<Dcl 5 (P) = Ai(P)« + A 2 

For Del <5 5 (P) = As(P) + A 4 (4) 

TABLE 1 

INSTRON 4206 COMPLIANCE CURVE FIT COEFFICIENTS 



Coefficient 


Value 


Units 


Ai 


.0512 


mm 


A 2 


-.0217 


mm 


n 


.1975 


— 


A 3 


.0064 


mm /kg 


A 4 


.0242 


mm 



B. BUNDLE TESTING 

The bundle testing was accomplished with the use of samples as 
illustrated in Figure 12. The copper tabs were prepared and washed 
with a dilute acid to enhance the performance of the adhesive. The 
graphite bundle was glued to the lower set of copper tabs with acety- 
lene adhesive. The upper tabs were then set so that the distance (d) 
between the the upper and lower tabs was the target gauge length. 
The graphite bundle was held straight with a 2 kg weight while the 
upper tabs were glued in place. The samples were carefully stored in 
the Composites Laboratory until use. As in the case of the compliance 
test sample, all of the samples tested were prepared by Laboratory 
Technician, Mr. Jim Nageotte. 
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Copper Tabs 




Figure 12 

Graphite Bundle Sample 



The bundle tests were performed on AS-4 Graphite bundles of 
3,000 filaments. The tests were run in accordance with the procedure 
outlined in Appendix D. The displacement due to compliance was sub- 
tracted from the total displacement and a summary of the the results 
is listed in Table 2. The data and test summary forms generated by the 
software are enclosed in Appendix E. 



35 



TABLE 2 



BUNDLE TESTING RESULTS (COMPLIANCE REMOVED) 



Material: 


AS-4 Graphite 


(3000 filaments) 






Test # 


Dry/ 


Gauge 


Modulus 


a 


P 


Comments 




Oil 




Length 




(mm/mm) 










(kg/mm) 








090901 


Diy 


5.0 cm 


1068 


3.214 


.0135 




090902 


Diy 


.5 cm 








* 


090903 


Dry 


.5 cm 








♦ 


090904 


Dry 


.5 cm 








♦ 


090905 


Dry 


.5 cm 








* 


090906 


Dry 


.5 cm 








★ 


090907 


Dry 


5.0 cm 








★ 


090908 


Dry 


5.0 cm 


912 


4.465 


.0145 




090910 


Dry 


2.5 cm 








4c 


080901 


Oil 


70.0 cm 






Test Validation 


100901 


Dry 


2.5 cm 


915 


3.000 


.0157 




100902 


Oil 


2.5 cm 








♦ 


100903 


Oil 


2.5 cm 








♦ 


100904 


Dry 


50.0 cm 


905 


3.035 


.0077 




100905 


Oil 


50.0 cm 


959 


3.660 


.0085 




100906 


Oil 


50.0 cm 








Not Analyzed 


100907 


Oil 


25.0 cm 


932 


3.964 


.0101 




100908 


Dry 


25.0 cm 


917 


3.180 


.0097 





* Adhesive Failure 



Examination of the load vs. percentage strain plots revealed the 
fact that for the shorter gauge lengths, 2.5 and 1 cm, the copper tabs 
failed to grip the graphite bundle uniformly. This caused some of the 
fibers to slip out of the tabs and the data from these tests were of no 
validity. The tests that had this phenomena are listed as adhesive fail- 
ure in the comments section of Table 2. Figure 13 shows a represen- 
tative example of this phenomenon as the load remains constant while 
the percent strain increases beyond the limit of the experiment. 
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Figure 13 




Load vs. Percent Strain for Sample 100903 (Adhesive Failure) 
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An examination of Table 2 shows that seven of the 18 samples 
resulted in satisfactory tests. The seven successful tests results have 
been summarized into Table 3. 



TABLE 3 

BUNDLE TESTING RESULTS (COMPLIANCE REMOVED) 

SUCCESSFUL TESTS 



Material; AS-4 Graphite (3000 filaments) 



Test # 


Dry/ 

Oil 


Gauge 


Modulus 

Length 

(kg/mm) 


a 


P 

(mm/mm) 


100901 


Dry 


2.5 cm 


915 


3.000 


.0157 


090901 


Dry 


5.0 cm 


1068 


3.214 


.0135 


090908 


Dry 


5.0 cm 


912 


4.465 


.0145* 


100908 


Dry 


25.0 cm 


917 


3.180 


.0097 


100907 


Oil 


25.0 cm 


932 


3.964 


.0101 


100904 


Dry 


50.0 cm 


905 


3.035 


.0077 


100905 


Oil 


50.0 cm 


959 


3.660 


.0085 



* Data not plotted 



Figures 14 through 32 are the graphs of load vs. percent strain of 
the test data, of the compliance, and of the test data with the fiber 
removed. 
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AS-+ BUNDLE TEST TEST DATA 



SAkPLUOOaOl GAJGE LENGTH » 25 MM 




Figure 14 

Load vs. Percent Strain Sample 100901 Test Data 
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Load CKq) 



AS-4 BUNDLE TEST COMPLIANCE 



SAMPLE 100901 C;(UGE LENGm = 15 MM 




Figure 15 

Load vs. Percent Strain Sample 100901 Compliance 
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PDO-I 



AS-4 BUNDLE TEST COMPLIANCE REMOVED 



SAliPLE 100^01 G4X;E length « 25 mm 




Figure 16 

Load vs. Percent Strain Sample 100901 Compliance Removed 
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Load CKaD 



AS-4 BUNDLE TEST TEST DATA 



SAIjPLE G9(}301 GA£E LENGTH » 50 UM 




Figure 17 

Load vs. Percent Strain Sample 090901 Test Data 
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Load CKaD 



AS-4 BUNDLE TEST COMPUANCE 



SAM>LE: 090301 GAJGE LENGTH » 50 MM 




Figure 18 

Load vs. Percent Strain Sample 090901 Compliance 
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Load C^^a^ 



AS-4 BUNDLE TEST COMPUANCE REMOVED 



SAtjPLC 090301 GfUGC LENGTH * £0 MM 




Figure 19 

Load vs. Percent Strain Sample 090901 Compliance Removed 
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AS-4 BUNDLE TEST TEST DATA 



SAMPLE 10030fi GAiyGE LEMJTH = 250 kW 




Figure 20 

Load vs. Percent Strain Sample 100908 Test Data 



45 



pdo-1 



AS-+ BUNDLE TEST COMPLIANCE 
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Figure 21 

Load vs. Percent Strain Sample 100908 Compliance 
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AS-4 BUNDLE TEST COMPLIANCE REMOVED 



SAtjPLE 10C905 GAUGE LENGTH » 2» MJ 




Figure 22 

Load vs. Percent Strain Sample 100908 Compliance Removed 
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Load C»^a!> 



AS-4 BUNDLE TEST (OIL) TEST DATA 



SWPLE 10M07 fiWlXJE LENGIH = 250 MM 




Figure 23 

Load vs. Percent Strain Sample 100907 Test Data 
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Load 



AS-4 BUNDLE TEST (OIL) COMPLIANCE 



SAMPLE 10M07 GAUGE LENGIH = 250 MM 




Figure 24 

Load vs. Percent Strain Sample 100907 Compliance 
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AS-+ BUNDLE TEST (OIL) COMPLNC REMOVED 



SAMPLE 10&307 GAUGE LOCTH = 250 MM 




Figure 25 

Load vs. Percent Strain Sample 100907 Compliance Removed 
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AS-+ BUNDLE TEST TEST DATA 



SMjPLC 1009(M GAiJGE LENGIH » 500 Myf 




Figure 26 

Load vs. Percent Strain Sample 100904 Test Data 
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Figure 27 

Load vs. Percent Strain Sample 100904 Compliance 
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AS-4 BUNDLE TEST COMPUANCE REMOVED 



SfllPLE 100904 GAIX3E LENGTH = 500 LW 




Figure 28 

Load vs. Percent Strain Sample 100904 Compliance Removed 
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Figure 29 

Load vs. Percent Strain Sample 100905 Test Data 
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Figure 30 

Load vs. Percent Strain Sample 100905 Compliance 
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Load 



AS-4 BUNDLE TEST (OIL) COMPLNC REMOVED 



SIMPLE 100^05 GAIUGE LENCTH 500 MM 




Figure 31 

Load vs. Percent Strain Sample 100905 Compliance Removed 
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V. CONCLUSIONS 



A. COMPLIANCE TESTING 

The curve-fitting procedure outlined in Appendix C was successful 
in arriving at a continuous function to describe the load-induced dis- 
placement of the load train for the INSTRON testing machine. Figure 
32 shows the data points measured during the zero gauge length test. 
Figure 33 shows the curve fit function using the parameters listed in 
Table 2. Figure 34 then shows the data points and the curve fit 
superimposed on one another. The greatest deviation in load is 
.0017 kg. 

Based on the small deviation of the data points to the curve fit 
function, the curve fit of the compliance is satisfactory to be used to 
subtract the displacement due to compliance from bundle test data. 

The curve fit solution for the machine compliance is satisfactory 
for loads greater than .1 kg. The curve fit does not do a good job in the 
region below this load. This is due to the inconsistent zero displace- 
ment length. The zero point for displacement is very difficult to 
achieve for bundle testing. The sensitivity of the load cell is insuffi- 
cient to consistently identify an actual zero load level on the bundle. 
The effect of this on the fiber strength distribution parameters is that 
some elements of the lower tail are distorted. 
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Figure 32 

0 Gauge Length INSTRON Compliance Test 
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Figure 33 

INSTRON Compliance Curve Fit Function 
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Figure 34 

0 Gauge Length INSTRON Compliance Test and Curve Fit 
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DISPLACEMENT (mm) (Dc - Do) 



The compliance testing method for the INSTRON has been 
demonstrated to be feasible and to provide realistic results. Testing 
should continue to determine whether the compliance is repeatable or 
whether it is a function of other factors such as temperature, humid- 
ity, or age of the machine. Full understanding of compliance is impor- 
tant because of the magnitude of compliance effects on short gauge 
lengths; it is the short gauge lengths to which stressing is applied. 

B. BUNDLE TESTING 

The bundle tests demonstrated the difficulty in obtaining test data 
for short gauge lengths. Table 2 illustrated the fact that no successful 
1 cm tests were completed and only one 2.5 cm test was completed. 
An improved technique of sample preparation to prevent the slipping 
of the fibers through the tabs is needed. Possibly an adhesive with bet- 
ter wetting qualities could be tried. 

The difficulty in testing the small gauge lengths highlights two 
important effects on prestressing. The first is that the shorter the 
prestress gauge length, the more difficult the process becomes due to 
the difficulty in maintaining gauge length. The second important effect 
is that the prestress effect on fiber parameter improvement will have 
to be based on a gauge length less than the ineffective length. The 
optimization process will now include another variable, the minimum 
possible gauge length due to machine constraints. As the minimum 
gauge length, increases the effectiveness of the prestress will 
decrease. 
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The importance of accurately removing the displacement in the 
test data for short gauge lengths was demonstrated by the bundle 
tests. Figures 15, 18, 21, 24, 27, and 30 show that as a function of 
percent strain, the shorter the gauge length, the greater the percent- 
age of test displacement is due to machine compliance. Figure 15 
shows that the percent strain displacement due to compliance for the 
shortest gauge length (2.5 cm) was greater than .6 percent, or nearly 
25 percent of the total test displacement. Figure 27 shows that the 
percent strain displacement due to compliance for the longest gauge 
length (50 cm) was less than .02 percent, which is less than 1 percent 
of the total test displacement. 

The inaccuracy incurred in not eliminating the compliance dis- 
placement for long samples before reducing the data to determine the 
fiber strength parameters is minimal. In the case of short samples, 
however, this is not the case. If the compliance displacement is not 
removed, the effect on the strength pcirameters will be the same as if 
the gauge length was actually longer than what was tested. This results 
in the fiber parameters being conservative in the prediction of fiber 
strength. While this is safe, it defeats the purpose of prestressing. 

Figures 25 and 31 are the graphs of the oiled bundles. Comparing 
these graphs to Figures 22 and 28, which are the dry or not oiled 
samples, the effects of friction can be seen in the shape of the curves 
after the maximum load has been reached. The dry sample load drops 
in large vertical jumps. These large jumps correspond to many fibers 



62 



breaking at once instead of sequentially. Figure 35 shows the 
simulated test data for a 100 mm sample. The curve is very smooth 
and does not show the sharp vertical drops the actual test data shows. 

The large vertical drops in the actual test data are noticeably 
reduced, but not eliminated, in the oiled samples. The oil. a soybean 
derivative, reduces the friction but does not eliminate all of the effects. 
The effect of friction on the test does appear to dominate the tail end 
or the final 30 percent of the percent strain. The effect of friction on 
the fiber strength parameters is difficult to predict. One effect that is 
clear is that the data for the first 50 percent to 60 percent of the 
breaks in the long samples do not appear to be affected by friction as 
the remaining breaks. In addition. Figures 16 and 19 indicate that the 
effect of friction declines as the gauge length becomes shorter. 

The need to eliminate friction between the fibers in a bundle test 
is important if long-gauge length tests are to be accurately completed. 
Different oils may be tried as lubricants to improve upon the perfor- 
mance of the soybean oil. If the friction cannot be eliminated, it must 
be brought down to a minimal level. 

The effect of noise in the bundle tests appears to be minimal. Fig- 
ures 36, 37, and 38 are graphs of simulated bundle tests with noise 
added to the load data. Figure 4 is a graph of the same simulation with 
no noise present. 

A noise level of 1 gm is equivalent to .002 percent for the 50 kg 
load cell used by the INSTRON. The difference in the curves from 
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Figure 35 

Bundle Test Simulation (Gauge Length— 10 cm) 
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Figure 36 

Bundle Test Simulation With ±1 gm Noise 
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Figure 37 

Bundle Test Simulation With ±10 gm Noise 
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BUNDLE SIMULATION DISCRETE DATA 

N = 1000 NOISE = +/- 10Q CMS 
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Figure 38 

Bundle Test Simulation With ±100 gm Noise 
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DISPLACEMENT (MM) 



Figure 35 to Figure 36 is barely noticeable. A noise level of 10 gms is 
equivalent to .02 percent. Figure 37 shows some effects of noise, but 
they are minimal. A noise level of 100 gms is equivalent to .2 percent. 
This is an excessive amount of noise and the effects of a noise level 
this high are clearly shown by Figure 8. 

The INSTRON manual specifies the sensitivity, or noise level, of 
the load cell to be .01 kg. This corresponds to roughly the noise level 
of Figure 37. An inspection of Figure 37 demonstrates that the effect 
of this noise level is barely discemable when compared to Figure 35, a 
simulation with no noise. 

This research has demonstrated that the physical process of 
preloading fibers with a short gauge length is difficult due to the 
physical limitations of the sample gripping mechanism. More work is 
required to find an adequate means to grip or restrain the fibers for 
short gauge lengths. 

The removal of the displacement due to compliance in the test 
results is significant for gauge lengths shorter than 2.5 cm in order to 
determine accurately the fiber strength distribution parameters. The 
effect of compliance diminishes with longer gauge lengths. The fric- 
tion between the fibers in bundle alters the shape of the load vs. dis- 
placement curve. Friction has more of an effect on longer samples 
than on short samples. Friction may be reduced by using oil on the 
samples, but it does not entirely eliminate the effect. Noise is not a 
significant factor in altering the data and does not interfere with the 
strength distribution parameters. 
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VI. RECOMMENDATIONS 



Further study into the magnitude of the machine compliance for 
the INSTRON 4206 testing machine should be initiated. The effects of 
temperature, the gripping devices, and repeatability should be investi- 
gated, To complete this, a quicker method to determine the curve fit 
coefficients needs to be developed. A Fortran program that can accu- 
rately and quickly solve the curve-fitting routine in Appendix C should 
be written and validated. A series of compliance tests may then be 
completed to build a database. 

The methods for sample preparation and bundle testing need to 
become less technician dependent. There are too many areas where 
mishandling of the samples may introduce errors into the data by 
damaging the fibers. The reliance on the laboratory technician to build 
the samples, align the samples, and then set the gauge length all add 
variability to the data that should be reduced if not eliminated. 
Hydraulic grips can eliminate the hand- and wrench-tightening pro- 
cedure required by the current grips. A brace or stand should be 
designed and built to hold the sample during loading into the 
INSTRON. This will eliminate much of the handling the sample must 
now endure. 

The prestress effects of preloading the fibers to a load greater 
than the critical load should be investigated. This is a probability 
problem whose solution could be integrated into the bundle simulation 
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program. TTiis program could then be used to predict the prestressing 
performance and then validated with tests performed in the laboratory 
on gauge lengths that can be successfully tested. This is very important 
because the poor test results of the short samples indicate that a short 
preload gauge length (< 1 cm) may be difficult to achieve. 

The effects of prestressing with respect to long-term reliability 
and composite fatigue failure should be investigated. The impact of 
local load stress cells on long-term reliability needs to be addressed 
quickly. The urgency of this aspect of prestressing is that due to the 
nature of fatigue failure, data taking is a very long process. 
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APPENDIX A 



MODEL SELECTION 

Composite fiber materials such as graphite, aramid, or boron are 
characteristically long and thin when compared to the dimensions of 
the composite. For example, the diameter of a graphite fiber is on the 
order of 5 microns, while the fiber length can extend to hundreds of 
feet for a 1 square foot pressure vessel. Under tension, the fibers fail, 
or break, at positions where flaws or imperfections exist. These flaws 
cause the fiber to break at loads much below the theoretical strength 
of the fiber, which is the chemical-bonding strength. The type of the 
flaw may be scratching, chemical erosion, contamination, or one of 
many others. The source of the flaw may have come in production, 
handling, or service use. The spatial location, density, and severity of 
these flaws determine the strength of each fiber. 

A single fiber or filament may be visualized as a chain made up of 
many links. These links, or segments, each have a unique strength 
that is limited by the severity of the flaws contained in that segment. 
The more severe the flaw, the weaiker the segment. For a long fiber, 
one with many segments, the strength of that fiber is only the 
strength of the most severe flaw, or the weakest link. If the severity of 
the flaw is known to fall within a certain probability distribution and 
the location of the flaws are randomly distributed along the fiber, then 
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a “weak link” model can be assumed for predicting the strength of a 
fiber. 

From statistics, the Weibull distribution best fits the “weak link” 
model. Phoenix and Wu, in Reference 2, describe the Weibull statisti- 
cal model as it applies to composite fiber strength prediction in detail. 
In general, the fiber is partitioned into a series of segments (m). The 
shortest limiting length segment is equal to the ineffective length or 
effective load transfer length of the fiber. This is on the order of mm 
for graphite/epoxy. A bundle is then formed by joining several, (n), 
fibers in parallel. For example, in this project, AS-4 bundles where m 
equals 1 and n equals 3,000 were tested in the laboratory. 

Based on the Weibull weak link model, the failure cumulative den- 
sity function (CDF) of a fiber can be explicitly written as a function of 
the two Weibull parameters a and p. The Weibull failure CDF is written 
as Equation 5. 



F = 1 - exp {-(x/P)a) (5) 

The two parameters, a and p, correspond to the distribution 
shape and scale, a adjusts the shape of the distribution by altering the 
skew and range of values. As shown in Figure 39, a high shape para- 
meter, or a, corresponds to a very narrow distribution; an a equal to 
3.5 approximates a normal distribution; and a smaller a corresponds 
to a very wide distribution with large tails. When a is large, the scale 
parameter p is approximately the mean. The AS-4 that was tested in 
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this project for a gauge length of 25 cm has an accepted a equal to 4 
and a |3 equal to 16 gm/ (mm/mm). 




Figure 39 

Effect of the Shape Parameter on the Weibull Distribution 

The Weibull distribution is not only a physically appropriate model 
for fiber strength, it has an additional advantage in that its failure 
cumulative distribution (CDF) can be written as an explicit equation. In 
addition, the CDF may be linearized with respect to |3. This results in 
Equations 6 and 7. These two equations may then be graphed together 
to form a graph such as the one shown as Figure 40. This graph of the 
F*/x* function allows direct graphical correlation of strength to failure 
probability. 
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Figure 40 

Weibull Linearized Failure Prediction Graph 

F* = ln(-ln(l - exp {-(P/p)"})) (6) 

X* = In (P/p) (7) 

where P is the load subjected to the filament. 

To predict reliability, which is one minus the failure probability 
for a fiber, the a and P must be known for that fiber to a sufficient 
degree of accuracy. Fiber bundle failure tests are performed in the lab- 
oratory to arrive at the parameters. A difficulty arrives with the fact 
that the parameters are a function of the length of a fiber. This corre- 
sponds to the fact that a long chain is weaker than a short chain 
because the longer chain has a higher probability of having a weak 
segment. 
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The Weibull model provides a means to normalize the parameters 
between different lengths of a common fiber if the parameters are 
known for any one length. This allows for the laboratory results to be 
applied to fibers of lengths different than the lengths tested. This fea- 
ture of the Weibull function is commonly called the size effect. 

If the strength distribution parameters pi and al are known for a 
fiber of length li and the parameters P2 and 02 are desired for the 
same fiber but of a different length I2 the following derivation provides 
the relationship: 



Fi = 1 - exp {-(x/pi)ai) 



( 5 ) 



R = 1 - F 



(Reliability = 1 - Failure) 



( 8 ) 



Ri = exp {-(x/Pi)“i} 



( 9 ) 



m = I1/I2 



( 10 ) 



Rt = n Ri R2 = (Ri)m 

Rt = [exp {-(x/pi)“i)p 



( 11 ) 



R2 = exp {-m(x/Pi)‘^i} 



( 12 ) 



R2 = exp {-(x/p2)“2} 
Equating Equations 9 and 12 then: 



( 9 ) 
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exp {-m(x/pi)“i} = exp {-(x/p 2 )“ 2 } 
{-m(x/pi)ai} = {-(x/P 2 )“ 2 } 
(x/(m)‘l/“l Pl)®l = (x/p 2 )^2 



( 13 ) 



Solving Equation 13 results in Equations 14 and 15. 



02 = ai 



(14) 



p2 = Pi(m)-i/«i 



(15) 



An implication of the size effect relations, Equations 14 and 15, is 
that a long composite made up of many segments may be reduced to a 
composite of one segment. This simplification provides for easier 
mathematical modeling of the composite as a large composite model 
shown in Figure 41 with many segments may be reduced to the much 
smaller model shown in Figure 42. This greatly reduces the amount of 
computational effort required to model composite strength or bundle 
strength as the model matrix has been reduced from a (m X n) to a (1 
Xn). 
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M = L/L* = 6 
N = 10 



Figure 41 

Composite Model Without Size Effect 
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Figure 42 

Composite Model With Size Effect 
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APPENDIX B 



BUNDLE FAILURE SIMULATION 

The purpose of the bundle failure simulation was twofold. The first 
was to allow the user to readily produce simulated bundle test data 
with variable fiber strength distribution parameters. The resulting data 
could then be graphically presented into load versus displacement or 
load verses percent strain representations. These graphs were utilized 
to gain understanding and insight into the effect of fiber strength 
parameters and gauge length on fiber failure process. 

The second purpose of the bundle failure simulation was to 
provide a means of verification for data interpretation and reduction 
software being developed for use in conjunction with the INSTRON 
Testing machine used in the Composites Laboratory, Naval 
Postgraduate School. The simulation could be used to ensure that the 
software returned the parameters from the simulated test data that 
were used as inputs to the bundle simulation. For this project, the 
compliance removal procedures and software described in Appendices 
C and D were validated with this program. Lt. Joseph Schmidt, in a 
separate project, used this program to validate software that reduced 
bundle test data into the fiber strength probability distribution 
parameters. 

The bundle failure simulation program “T2” takes fiber strength 
distribution parameters a and P for the length of the bundle, the 
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number of filaments, the length of the bundle, the fiber modulus, and 
testing machine settings of cross-head speed and sampling rate, and 
then produces two sets of load vs. displacement data. The first set is 
the loaded and displacement determined for each individual filament 
break. This is referred to as the analog data. The second set, which 
simulates fiber testing machine data, is the continuous load and 
incremental displacement data, which is referred to as discrete data. 
The program also has the added feature of inputting a noise range. The 
maximum number of filaments is 1,000. 

The program solves for the strength of each filament based on the 
Weibull distribution and the size effect. An explanation of the Weibull 
function and of the size effect is detailed in Appendix A. 

A. ANALOG DATA SIMULATION 

The analog data was determined by using a random number 
generator to to generate a random strength of each fiber, determine 
the displacement of each break, and then compute the load of the 
bundle at each break. 

Equation 16 was used to determine the strength of each fiber. 

Pf = ln(-ln(l - exp {-(x/|3)a})) (16) 

where Pf = failure load of the the fiber 
X = random number 0 < x < 1 

Equation 17, which is a manipulation of Hooke’s law, is used to 
determine the displacement of each fiber at the failure load. 
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df = Pf * 1/E 



(17) 



where df = Displacement of the Fiber 

1 = Length of the Fiber (Gauge Length) 

E = Strain Modulus 

The displacement at each break was stored in an array and then 
sorted in ascending order. The load of the bundle at each break was 
then determined by using Equation 18. This is also a manipulation of 
Hooke’s Law. At each break, because the simulation is for a controlled 
displacement, the load in the bundle drops by a factor of 1/n. This 
data point is solved for by Equation 19. 



where i = 1,2,3, . . . n (Number of Breaks) 

B. DISCRETE DATA SIMULATION 

The discrete data was found by utilizing the analog data. The 
analog data was transformed into a continuous function by connecting 
the data points. The displacement between data points is found by 
Equation 20. 



(Pb)i = ((df)i * E * (n-i +1))/1 



(18) 



(Pb)i+ = ((df)i * E * (n-i))/l 



(19) 



Deld = CHS/SR 



( 20 ) 
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where CHS = Cross-head Speed (mm/min) 

SR = Sampling Rate (pts/sec) 

The discrete data displacement is in increments of Deld while the 
discrete data load is the intercept of the analog continuous data for 
that appropriate displacement. Equation 21 is used to compute the 
discrete displacement points. Discrete data points are solved for until 
the load reaches zero. 

d = Deld * 1 1 = 0,1.2, . . . Load = 0 (21) 

C. NOISE SIMULATION 

The user has the option of simulating the effects of noise in the 
load transducer of a bundle test. If desired, the noise is added to the 
discrete load data by Equation 22, 

Pnolse = P - ((Ph - Pl)/2) + X * (Ph - Pi) (22) 

where P = Original Load 

Ph = High Noise Range 
Pi = Low Noise Range 

D. PROGRAM LISTING 

The program was written in Microsoft Fortran 4.01. It is designed 
to be operated from any IBM/ Compatible computer with two floppy 
drives or hard disk and one floppy drive. The executable program is 
run from the A; drive or the hard disk and the data is written to the B: 
drive. To execute the program type: T5 <enter>. 
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SUBROUTINE LISTING: 



INIT 


Initializes arrays to zero 


INPUT 


Prompts user to input required value 


RAND 


Random number generator (user must input seed) 


SLOAD 


Solves for the breaking stress of each filament 


SORT 


Sorts the filaments by breaking stress and then 
solves for breaking displacement 


ANLDAT 


Creates the analog load and displacement file 


DISDAT 


Solves for the continuous and creates the discrete 
load and displacement file 


NOISE 


Allows user to input random noise to a specified 
threshold into the discrete load data 


OUTPUT 


Writes the analog and discrete data to floppy 
disks 



VARIABLE LISTING: 

DATA(0: 1000,7) Workspace array 



ADATA(0:2000,2) 


Analog data output file; Displacement, Load, listed 
sequentially by first break to last 


DDATA(0:2000,2) 


Discrete data output file; Displacement, Load, 
listed as a function of sampling time 


ALPHA 


Weibull shape parameter for the given length 


BETA 


Weibull scale parameter for the given length 


E 


Modulus of fiber 


L 


Length of the Fiber 


N 


Number of filaments in the bundle (max. = 1,000) 


CHS 


Simulated testing machine cross-head speed 
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SR 


Simulated machine sampling rate 


DELD 


Distance between data points for discrete data 


POINTS 


Number of discrete data points 


SEED 


Random number generator seed for strength 
computations 


SEEDl 


Random number generator seed for discrete data 
load during break generator 


SEED2 


Random number generator seed for noise 
generator 


BREAKS 


Counter that is used to find slope of analog data 
for up to n- 1 breaks 


LIMIT 


Maximum strength of fiber i used in solving for 
the maximum displacement for n-1 breaks 


PH 


Upper noise range 


PL 


Lower noise range 


OUTPUT FILES 




ADATA 


Analog data; displacement, load 


DDATA 


Discrete data; index, displacement, load 
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PROGRAM T2 



* PROGRAM T2 WILL SIMULATE A BUNDLE TEST FOR N FIBERS IN PARALEL 

* FOR GIVEN FIBER STRENGTH PARAMETERS AND LENGTH OF THE BUNDLE. 

* THE PROGRAM WILL OUTPUT ANALOG DATA THAT IS A LISTING OF THE * 

* DISPLACEMENT AiND LOAD OF EACH FIBER BREAK. THE PROGRAM WILL 

* OUTPUT DISCRETE DATA THAT SIMULATES MACHINE TESTING DATA FOR THE * 

* BUNDLE. THE OUTPUT IS BY INDEX, DISPLACEMENT AND LOAD. THE 

* NUMBER OF DATA POINTS IS A FUNCTION OF THE SIMULATED CROSSHEAD 

* SPEED AND THE SAMPLING RATE. * 

* Vf 



* THE OUTPUT IS WRITTEN TO DRIVE B: ADATA. AND DDATA. * 

•>v * 



Vf * y? 7 V * ?V * "s V Vr Vr ' 



*5 V Vc * - jV 7’c Vc Vr * V? Vr >V 7V V? 7V 7 V 7V 7V • 



’ 7V >V 7V 7V 7V 7V 7V 7V 7V 7V 7V 7V 7V 7V 7V 7V 7V Vc 7V 7V 7V 



’ 7V 7V 7V 7V 7V 7V 7V 7V 7V Vc 7 V 7V 



* DECLARATION STATEMENTS 

REAL DATA(0; 1000,7), DDATA(0; 5000,2) , ALPHA, BETA, L, CHS, PH, PL, 
1E,DELD,SR,ADATA(0; 2000,2) 

INTEGER N, I, J,K, POINTS, Q 
1 CONTINUE 



* PRINT BANNER AND START PROGRAM 

PRINT*, 'COMPOSITE BUNDLE TEST SIMULATION PROGRAM’ 

CALL INITCDATA, DDATA, ADATA) 

CALL INPUT( ALPHA , BETA , L, CHS . E , N , SR ) 

CALL RAND(DATA,N) 

CALL SLOAD( DATA. N, ALPHA, BETA) 

CALL SORT(DATA,N,E,L) 

CALL ANLDAT(DATA,N) 

CALL D I SDAT ( DDATA , DATA ,N,E,L, CHS, SR, POINTS) 

CALL NOISE(DDATA,POI.NTS) 

CALL OUTPUTC DATA , DDATA , N , E , L , ALPHA , BETA , POINTS , ADATA ) 

* ALLOW USER TO RERUN PROGRAM 

PRINT* 

PRINT*, 'RUN COMPLETE, DO YOU WISH TO RUN ANOTHER SET OF DATA' 
PRINT*,' IF YES TYPE 1, ELSE TYPE 0 TO EXIT PROGRAM’ 

READ*, Q 

IF(Q .EQ. 1) GOTO 1 

* PRINT BA.NNER AND END PROGRAM 

PRINT* 

PRINT*,' END OF RUN, HAVE A NICE DAY' 

END 

* SUBROUTINE INIT 

* SUBROUTINE INIT INITIALIZES THE ARRAYS TO ZERO 



* DECLARATION STATEMENTS 

SUBROUTINE INIT( DATA , DDATA , ADATA) 

REAL DATA( 0; 1000 , 7 ) , DDATA( 0; 5000 , 2 ) , ADATA( 0; 2000 , 2 ) 
INTEGER I,J 

* PRINT BANNER AND START SUBROUTINE 
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